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A  linear  structure  can  be  characterized  by  its  compliance  matrix 
C,  which  is  6x6  symmetrical  and  positive  definite  and  which  relates  a 
force  6-vector  F  .  to  a  displacement  6-vector  D.  by  the  relation  CF  .  ■  D.. 

j         .  J  j      j 

The  inverse,  S,  of  C,  is  called  the  stiffness  matrix  and  satisfies  SD .  ■ 

J 

P  .   This  thesis  deals  with  the  problem  of  finding  optimal  values  of  such 

matrices  C  and  S  from  experimental  determinations  of  a  sufficient  number 

of  vector-pairs  (F..D.)  which  are  presumed  to  contain  random  errors. 

J   J 

J.  E.  Brock  has  introduced  this  problem  area,  suggested  several  dif- 
ferent criteria  of  optimality,  and  solved  some  of  the  corresponding 
specific  problems.   This  thesis  completes  the  solution  to  a  previously 
unsolved  specific  problem  of  this  group  and  contributes  computationally 
convenient  new  solutions  to  another.   Moreover,  a  computer  program, 
originally  written  for  the  CDC  1604  has  been  rewritten,  in  FORTRAN  IV 
Language,  as  two  programs  for  the  IBM  System  360  computer,  and  the 
capability  has  been  significantly  augmented. 
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SECTION  1 
INTRODUCTION 
1.1  Source  and  Statement  of  the  Physical  Problem. 

The  coordinate  system  used  here  is  a  right  handed  orthogonal  triad, 
whose  unit  vectors  are  denoted  by  e. ,  e  ,  and  e  ;  see  figure  1. 
A  generalized  force  will  be  denoted  by  F,  where 

F«  {Fr  F2,  F3,  F4,  F5,F6  } 

is  a  column  matrix.   Braces  £}■  indicate  that  a  column  matrix  is  in- 
tended.  The  first  three  elements,  F- ,  F_,  and  F  ,  are  force  components 
in  the  e. ,  e  ,  and  e  direction  respectively,  and  the  last  three  ele- 
ments, F,,  F  ,  and  F,,  are  moment  components  about  these  axes. 
Also  let  a  generalized  deflection  be  denoted  by  D,  where 


B  =  iDi>  WW  D6] 


The  first  three  elements,  D1 ,  D«,  and  D~,  are  displacements  in  the 
e. ,  e  ,  and  e_  direction  respectively,  and  the  last  three  elements  are 
rotation  components  about  these  axes. 

It  is  now  presumed  that  a  force  F  is  applied  to  a  (suitably  con- 
strained) elastic  structure,  the  application  being  at  point  P,  and  the 
corresponding  deflection  D  is  observed  at  P;  see  figure  1. 

It  is  postulated  that  a  linear  relation  holds  between  the  elements 
of  F  and  D,  i.e.,  that  there  exists  a  matrix  C  such  that 

CF  -  D 
It  is  further  postulated  that  matrix  C  is  symmetric  and  positive  definite, 
and  that  it  has  an  inverse  matrix  S  (S  =  C  )  which,  likewise,  is  symmetric 
and  positive  definite. 


FIGURE  1 


Matrix  C  is  called  the  compliance 
(or  flexibility)  matrix  (of  the  structure) 
and  its  elements  are  called  compliance  or 
flexibility  coefficients. 

Matrix  S  is  called  the  stiffness  matrix, 
and  its  elements  are  called  stiffness  in- 
fluence coefficients. 

Note  that  matrices  C  and  S  depend 
not  only  on  the  structure,  but  also  on 
the  choice  of  point  P. 

There  are  basically  two  approaches  for  obtaining  the  compliance 
matrix  C,  or  the  stiffness  matrix  S.   One  employs  elastic  theory  in  an 
attempt  to  come  up  with  the  elements  of  C  or  S.   Even  with  simple  struc- 
tures it  may  not  be  easy  to  write  down  the  elements  of  these  matrices. 
The  theory  must  be  stretched  beyond  certainty  and  one  obtains  only  an 
approximation  to  the  correct  result. 

The  second  approach,  the  use  of  which  leads  to  this  thesis  problem, 
employs  experimental  data  rather  than  elastic  theory  to  come  out  with  the 
best  C  or  S.   Because  of  inevitable  random  errors  present  in  the  data,  it 
is  impossible  to  obtain  the  matric  C  or  S,  without  using  some  mathemati- 
cal procedure  that  is  able  to  minimize  the  influence  of  these  errors, 
and  at  the  same  time  to  assure  that  matrices  C  and  S  are  symmetrical  and 
positive  definite.   Furthermore  the  data  must  contain  a  sufficient  num- 
ber of  independent  pairs  of  vectors  F  and  D. 

The  question  of  how  many  pairs  of  vectors  are  required  is  consider- 
ed in  section  7.   However  in  practice,  it  is  anticipated  that  many  more 


pairs  will  be  available  than  the  theoretical  minimum. 


1.2  Brief  History  of  this  Problem  and  Related  Problems. 

Professor  John  E.  Brock  [1],   of  the  Mechanical  Engineering  Depart- 
ment at  the  Naval  Postgraduate  School,  introduced  this  problem.   He  has 
defined  different  criteria  of  optimality  and,  corresponding  to  each  such 
criterion,  has  established  a  matrix  equation,  the  solution  of  which  will 
minimize  the  corresponding  errors  and  permit  determining  the  desired  sym- 
metric and  positive  definite  C  or  S  matrix. 

In  the  next  chapter  there  will  be  given  the  mathematical  descriptions 
of  four  criteria  of  optimality  that  Professor  Brock  has  introduced  in 
reference  [1],  and  the  developments  leading  to  the  matrix  equations  that 
which  must  be  solved  for  an  unknown  square  matrix. 

A  description  and  an  explanation  of  how  to  use  the  computer  programs 
which  has  been  developed  for  solving  these  equations  will  be  given  later. 

The  data  used  to  test  the  programs  are  not  from  an  actual  experi- 
ment.  An  explanation  will  be  given  of  how  sample  data,  containing  random 
errors,  are  generated  for  test  purposes. 

In  this  thesis,  matrix  C  has  been  restricted  to  be  symmetrical  and 
positive  definite,  as  is  required  in  the  structural  problem  under  con- 
sideration.  Two  other  related  situations  will  be  considered  in  Appendix 
D. 

In  general  the  different  criteria  of  optimality  that  will  be  given 
in  the  next  chapter  can  be  applied  to  all  systems  which  are  described  by 
a  matrix  equation  of  the  form  CF  =  D,  where  C,  F,  and  D,  need  only  to 
be  comformable  for  matrix  operations. 


Numbers  in  brackets  refer  to  references  listed  in  the  Bibliography 
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SECTION  2 
MATHEMATICAL  DEVELOPMENT  OF  THE  PROBLEM 
2.1  Physical  Statement  of  Criterion  One. 

It  is  assumed  that  there  are  n  pairs  of  vector  observations  F  and 
D. 

As  criterion  one,  it  is  desired  to  determine  C  so  as  to  minimize 
the  sum  of  the  square  of  all  scalar  errors,  viz.,  to  minimize 
-n  6  i 


where  {  }  .  =  L-th     element  of  vector  | J 


(1) 


Now  let 


0  -   [  ?v    F2,  F3,  ,  Fn] 


V  =   [  Dr,  D2,  D3,  ,  Dn] 


These  are  (  £  *  'T?   )  matrices. 
From  equation  (1) 


o<,--|rc{^-5J.}>c{cFJ-5J.}2f+...+(^.-5i]j]    (2) 


Expanding  the  first  term  of  equation  (2) 


I  (i#rW5 = C{^.-*W  +  (Icp.-U)  +.:.  +  a^-M,  )        o) 


i8' 


If  all  the  terms  in  equation  (2),  are  expanded  as  in  equation  (3), 

it  is  easy  to  see  that  these  terms  are  nothing  other  than  the  diagonal 

T 
terms  of  the  matrix  product:  (CU-V)  (CU-V) ,  and  therefore  it  is  possible 

to  write 
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OC,  =  tn.   J3(CU-V)T(CU-V)] 


(4) 


where  tr  stands  for  the  trace  of  a  matrix  and  superscript  T  denotes  the 
transposed  of  a  matrix. 

2.2  Alternate  Criteria. 


J>    A 


>  F 


FIGURE  2 
As  may  be  seen  from  Figure  2,  in  a  one  dimensional  case,  criterion 
one  requires  minimizing  the  sum  of  the  square  of  the  vertical  errors  a.. 
It  is  natural,  therefore  to  consider,  as  criterion  two,  minimizing  the 
sum  of  the  squares  of  the  horizontal  errors,  b..   In  the  six  dimensional 
case  that  this  thesis  is  concerned  with,  this  means  minimizing 


(5) 


Equation  (5)  is  in  the  same  form  as  equation  (1),  with  the  roles  of 

F,  and  D.  interchanged  and  using  the  stiffness  matrix  S  instead  of  its 

J      J 
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inverse,  the  compliance  matrix  C. 

Using  the  same  approach  as  before  this  sum  can  be  written  in  the 
following  way: 


<^%^iji  [(sv-u)T(sv-u)] 


(6) 


It  is  possible  too,  to  minimize  the  sum  of  the  areas  of  the  "error 
triangles",  whose  legs,  in  a  one  dimensional  case,  are  a.  and  b  ;  this 
is  the  same  as  minimizing 

<X3  =  tn   [(CU-V)T(U-SV)]  (7) 

Now,  think  also,  of  minimizing  a  weighted  sum  of  the  square  of  the 
horizontal  and  vertical  errors  (  /3,  °( ,  +y3r°(o  )•   This  criterion 
number  four  is  the  same  as  minimizing 

oCH   -t^[ACcU-v/(CU-V;H-y93.(SV-U)T(5V-U)]    (8) 

where  £  %   and  /3«  are  scalar  quantities  which  may  be  specified  so  as  to 
designate  the  relative  importance  of  the  two  separate  parts  of  which  it 
is  composed. 

From  equation  (8),  without  loss  of  generality,  it  is  possible  to 
make   /5 ,   =  1  and  J2>  2  =  Q>         ,  and  write  equation  (8)  as 

<*h  =1>l[(cu-V)T(cU-V)+^3(SV-U)T(SV-U)]     (9) 

Brock  [1]  also  suggested  a  fifth  criterion  of  optimality  but  since 
no  progress  has  been  made  toward  satisfying  this  criterion,  it  is  not 
included  here. 

2.3  Five  Theorems  on  the  Trace  of  a  Matrix. 

In  deriving  the  matrix  equation  corresponding  to  each  criterion  of 
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optimality,  the  following  five  simple  theorems  concerning  the  trace  of 
a  matrix  will  be  used.   The  proofs  are  simple  and  are  omitted. 

THEOREM  1:   If  the  matrix  P  and  the  matrix  Q  are  conformable  in 
either  order,  then  each  product  is  square  and 
tr  (PQ)  =  tr  (QP) 
THEOREM  2:   If  the  square  matrices  P  and  Q  may  be  added,  then 

tr  (P)  +  tr  (Q)  =  tr  (P  +  Q) 
THEOREM  3:   If  P  is  a  square  matrix,  then 

tr  (P)  -  tr  (PT) 
THEOREM  4:   If  P  is  a  matrix  and  x  is  a  scalar,  then 


THEOREM  5:   If  P  is  a  square  matrix  and  x  is  a  scalar,  then 

Note  that  theorem  1  permits  one  to  write: 

tr  (ABCD)  =  tr  (BCDA)  =  tr  (CDAB) ,  and  so  on. 


2.4  Derivation  of  the  Matrix  Equations  which  assure  the  Optimization. 

In  sections  2.1  and  2.2  four  different  criteria  of  optimality  have 
been  given.   This  section  devotes  attention  to  the  derivation  of  a  matrix 
equation  corresponding  to  each  criterion.   Solving  any  such  equation 
assures  achieving  the  corresponding  optimization. 
a)  Equation  (4)  is 

a,  =&.[(cu-v)(cu-v)] 

=&.(uTccu-vTcu-uTcv+vTv) 

T 
Note  that  instead  of  C  ,  C  has  been  written  since  C  is  symmetric. 

-\  r 

Now  let  x  be  any  given  element  of  C,  and  let  C*  denote   "  ~  ,  then 
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1*L     - 
3X 


=  fo(uTc*cu+uTcc*u-vTc  u-uTc*v) 


r&i[(cUUT^UUTC-UVT-VUTj(C*)] 

**i[(cp  +?c  -q)(c*)] 


where 


P  =  UUT,  and  Q  =  UV  +  VUT 


are  obviously  symmetric. 

Thus,  if  C  satisfies  the  matrix  equation 

CP  +  PC  =  Q  (10) 

then    -r — ' —  ■  0,  for  all  elements  x  of  C  and  there  is  an  extremum, 
which  will  be  minimum,  since  oC  |  is  tne  sum  °f  tne  squares  of  real 
numbers . 

Because  of  the  simple  relationship  between  oCt   and  o(  ^  ,  it  is 
evident  that  a  similar  equation  applies  for  criterion  two. 

This  equation  is 


where  0  is  as  before  and  R  =  W 
b)   Equation  (7)  is 


SR  +  RS  =  Q  (11) 

T 


o<3  -in  [(CU-V)  (U-SV)] 


The  solution  will  be  to  determine  the  matrices  C  which  make  o( * 
stationary;  among  these  will  be  one  which  minimizes   o(  3 • 

Knowing  that  S  =  C   ,  and  both  C  and  S  are  symmetric,  after  some 
manipulation 
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and  S*  = 


o(3  -U.  (UTCU-VTU-UTV+Vn5v) 

>efore,  let  x  be  any  given  element  of  C,  and  let  C*  =  4 — 


Since  CS  =  I,  then 
CS*  +  C*S  =  0 
and        S*  =  -SC*S. 
Thus 


|fi  =  u  (uTc*u  -vxs  c*s  v) 
=  k[(uuT-svvTs)(c*)] 


This  will  vanish  and  oC^  will  be  extremal  if  S  satisfies  the  matrix 

equation 

T  T 

UU     =  SWS 

After  post  and  pre-multiplying  by  C 

CUUTC  =  WT  (12) 

c)   Equation  (9)  is 

°CH  =  fa[(cu-v)T(cu-v)  +j3  (sv-uJT(sv-u)J 
*&[vTv-uTcv-vTcu  +UTCCU  + 
yS(vTss  v  -  Jsv  -  vTsu  +  uTu)] 

As  before,  let  x  be  any  given  element  of  C. 

Thus,  using  notation  introduced  previously,  viz.  : 


P  =  UUT,   R  =  WT,  and  Q  =  VUT  +  UVT 


it  may  be  seen  that 
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§|^fc[-uTrv-vc  u  +  uTc*cu+uTcc*u  + 

p  (vVsv  +vTss*v  -  uTs*v  -VT5*U)] 
=  fe[(-vuT-uvT+cuuT+uuTc)(c*)  + 
y3  (5vvT+vvTs  -  V(JT-  uvT)(s*)] 

=  ti  [(-a  +  cp+pc)(c^v3(SR+RS-fi)f5C^] 

=  tn  [  {  (-  4+PC+CP)  -J3S  (sR+RS-Q)  s  }  (c*)  ] 

which  will  vanish  if 

PC  +  CP  -  Q  =^3  S(SR  +  RS  -  Q)S       (13) 
If  the  right  hand  side  of  equation  (13)  is  set  equal  to  zero,  the 
left  hand  side  is  precisely  the  same  as  equation  (10) ,  which  corresponds 
to  criterion  number  one.   If  the  left  hand  side  is  set  equal  to  zero, 
the  right  hand  side  is  nothing  other  than  equation  (11),  which  corres- 
ponds to  criterion  number  two. 
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SECTION  3 
SOLUTION  OF  EQUATIONS  FOR  CRITERIA  1  AND  2 

Equations  (10)  and  (11)  that  minimize  criteria  1  and  2  respectively 
are  of  the  same  form,  where  the  matrices  C,  S,  P,  R,  and  Q,  are  all  sym- 
metric.  The  solution  by  Brock  [1],  given  here,  is  for  equation  (10),  but 
is  applicable  to  equation  (11),  with  the  roles  of  P  and  R  interchanged. 
The  solution  in  equation  (10)  will  be  for  the  unknown  matrix  C,  where 
in  equation  (11)  is  for  the  unknown  matrix  S. 

In  the  physical  problem  of  primary  interest,  matrices  C,  P,  and  Q 
are  all  6x6,  but,  for  the  sake  of  generality,  it  will  be  presumed  that 
they  are  m  x  m. 

In  equation  (10) 

CP  +  PC  =  Q 
there  are  really  only  N  =  %  m  (m  +  1)  unknown  elements  of  C.   These  un- 
knowns can  be  arranged  in  a  column  matrix  according  to  the  following 
procedure: 

Use  the  notation  M.  =  %  (j  -  1)  j. 

Given  i,  j,  let  t    =  Mmav/<  4\  +  min  (i,j). 
i,j    max(i,j; 

For  example,  with  i  =  3,  and  j  =  5,  then 

C3,5  =  C5,3  =  %  (4)(5)  +  3  =  13 
and  t  fi  =  t    -  h   (5)(6)  +  2  =  17 

Table  1  gives  values  of  t.  .  for  i,j  =  1,  2,  3,..., 10 

Next, a  column  matrix  W  of  N  elements  is  formed  from  the  N  unknown 

elements  c. .  of  matrix  C.   Element  c. .  of  C  becomes  the  t.  .th  element 

of  W;  note  that  t    =  N. 

m,m 


17 


i-1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

J-l 

1 

2 

4 

7 

11 

16 

22 

29 

37 

46 

2 

2 

3 

5 

8 

12 

17 

23 

30 

38 

47 

3 

4 

5 

6 

9 

13 

18 

24 

31 

39 

48 

4 

7 

8 

9 

10 

14 

19 

25 

32 

40 

49 

5 

11 

12 

13 

14 

15 

20 

26 

33 

41 

50 

6 

16 

17 

18 

19 

20 

21 

27 

34 

42 

51 

7 

22 

23 

24 

25 

26 

27 

28 

35 

43 

52 

8 

29 

30 

31 

32 

33 

34 

35 

36 

44 

53 

9 

37 

38 

39 

40 

41 

42 

43 

44 

45 

54 

10 

46 

47 

48 

49 

50 

51 

52 

53 

54 

55 

TABLE  1.   Values  of  t 


i.J 


The  column  matrix  Z  of  N  elements  Is  formed  in  a  similar  way  from 
the  known  elements  q.  .  of  symmetric  matrix  Q. 


where  z 


and   d    is  the  Kronecker's  delta. 

Finally  the  coefficient  matrix  P  is  formed  as  the  sum  of  m  (N  x  N) 

matrices  P_,  I  ■  l,2,...,m,  each  of  order  N  x  N,  which  themselves  are 

T 
formed  from  the  (m  x  m)  matrix  P  ■  UU  ,  according  to  the  following  pro- 
cedure: 

P  is  formed  from  P  by  inserting  rows  and  columns  of  zeros,  an  oper- 
ation that  is  denoted  by  (*)(K)(L)  which  means  insert  L  rows  of  zeros 

after  the  K   row  in  the  original  matrix  P,  and  insert  L  columns  of 

feh 
zeros  after  the  K   column  in  the  original  matrix  P.   The  P  is  obtained 

from  P  by  the  operation: 
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(a) 

waooij) 

(b) 

(*)(i)(i-d 

(c) 

(*)(i+i)(D 

(d) 

(*) (1+2) (I+l) ,   etc 

until  such  insertions  fall  outside  the  range  of  the  N  x  N  size. 

Then: 

m 

1=1 
and  solving  the  system 

PW  =  Z 

gives  the  column  matrix  W,  whose  elements  can  then  be  rearranged  to 

produce  the  desired  matrix  C. 

Bickley  and  McNamee  [4]  have  considered  a  more  general  equation  of 
the  form: 

AZ  +  ZB  =  F 
and  they  give  four  methods  of  solution,  which  they  have  named:  The  Big 
Matrix,  the  Irrational  Solution,  the  Semi-Rational  Solution,  and  A  Ration- 
al Solution. 

The  method  given  by  Larsen  [5],  and  the  method  given  by  Brock  [1] 
and  reproduced  here,  appear  to  be  additional  distinct  methods;  both, 
however,  probably  are  closely  related  to  what  Bickley  and  McNamee  refer 
to  as  the  "big  matrix",  in  which  A  ■  B. 

As  the  results  obtained  with  the  method  given  here  is  quite  satis- 
factory, no  further  study  has  been  done  with  the  methods  given  in  refer- 
ence [4],  but  it  might  be  of  interest  to  one  working  with  this  problem 
to  try  suitably  specializing  Bickley  and  McNamee's  methods. 
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SECTION  4 

SOLUTION  OF  EQUATION  FOR  CRITERION  3 

T 
To  solve  equation  (12),  postmultiply  by  UU  ;  then 

CUUTCUUT  =  WTUUT  (14) 

CUUT  =  (WTUUT)^  (15) 

T       T 
In  equation  (14),  UU  and  W  ,  both  are  symmetric,  but  the  product 

T  T 
W  UU   is  non-symmetric. 

Brock  [6],  has  solved  equation  (15),  using  a  highly  sophisticated 

computer  library  subroutine,  based  on  a  Laguerre-Hessenberg  algorithm, 

T  T 
to  find  eigenvalues  of  W  UU  ,  which  are  then  used  in  a  method  described 

by  Frazer,  Duncan,  and  Collar  [3],  in  order  to  find  CUU  =  (W  UU  )  , 

from  which  C  may  easily  be  obtained.   This  procedure  is  available  in  a 

computer  program  for  the  FORTRAN  63  language  to  be  used  in  the  CDC  1604 

computer. 

Owing  to  the  difficulties  in  converting  this  program  to  FORTRAN  IV 
language  used  in  the  IBM  System  360  computer,  recently  installed  at  the 
Naval  Postgraduate  School>  the  following  two  other  procedures  were  develop- 
ed and  have  successfully  been  applied  to  solve  equation  (12). 

a)   Iterative  procedure  for  finding  the  square  root  of  a  non- 
symmetric  matrix. 

Pulay  [7],  has  given  a  method  that  determines  the  principal  square 
root  of  a  positive  definite  symmetric  matrix  using  an  iteration  pro- 
cedure. 

Using  his  idea  with  some  modification,  a  method  for  obtaining  the 
square  root  of  a  positive  definite  non-symmetric  matrix  has  been  suc- 
cessfully used  and  applied  to  several  examples. 

In  equation  (15) 
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T      T   T  Jc 

CUU  =  (WUUT 


T  T 

let  WUU  =  A 


and  A^  -  D  +  B  (16) 

where  D  is  a  diagonal  matrix,  whose  elements  are  the  positive  square 
roots  of  the  diagonal  elements  of  A,  and  B  is  a  matrix  formed  from  CUU 
D,  and  C  is  an  approximation  of  the  matrix  C  that  will  minimize  the  cri- 
terion number  3. 

ft* 

The  matrix  C  can  be  obtained  from  the  equation 

CU  =  V 

Ifl  rp 

CUU    =  vu 

rS  T         T    -1 

and  C  =  VU   (UU   ) 

This  is  one  way  of  finding  the  matrix  C,  suggested  by  analogy  to 
Pulay's  method.  However,  in  the  computer  program  used  in  this  thesis 
to  solve  equation  (15),  the  matrix  2  is  taken  to  be  the  matrix  C  that 
satisfies  criterion  number  2.  This  is  done  in  order  to  save  computer 
time. 

From  equation  (16) 

(D  +  B)2  =  A 

DB  +  BD  =  A  -  D2  -  B2 

2 
This  equation  can  not  be  solved  for  matrix  B  due  to  the  term  B  , 

but  following  Pulay's  procedure  it  can  be  transformed  into  a  simple 

iteration  formula 

DBk  +  BkD  =  A  -  D2  -  B2^ 

and  the  (i,j)-th  element  of  B  ,  can  be  obtained  in  a  simple  manner  from 


the  formula 


(^■--^[(^VOUJ 
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To  obtain  (R  ) . .,  the  iteration  procedure  must  be  carried  out  until 
the  difference  between  two  successive  approximations  of  (B  ) .  .  become 
less  than  a  specified  allowable  error,  i.e. 

(B.).  .  -  (B,  .)..=  A 
k  ij     k-1  ij 

where  ^  in  the  actual  program  is  10   ,  and  the  ()  .  implies  comparing 

element  by  element. 

After  matrix  B  has  been  obtained 

C  =  (D  +  B)(UUT)_1 
or  C  =  (WTUUT)^(UUT)'1 

The  procedure  described  above  leads  to  output  3A  in  the  actual 

computer  program.   Output  3B  is  the  solution  of  basically  the  same  equa- 

T 
tion,  but  instead  of  equation  (12)  being  postmultiplied  by  UU  ,  it  is 


premultiplied  resulting  in 


T        T   T  h 

UU  C  =  (UU  w  y 


T  T 
and  using  the  same  procedure  as  before  the  square  root  of  UU  W  can  be 

obtained,  and  therefore  the  unknown  matric  C. 

Output  3C  is  simply  the  arithmetic  average  of  3A  and  3B. 

b)  The  second  procedure  that  solves  equation  (12)  uses  some 

matrix  algebra  manipulation,  and  instead  of  working  with  a  non-symmetric 

matrix  ,  it  requires  taking  the  square  roots  of  two  different  symmetric 

matrices.   This  procedure  was  called  to  the  attention  of  the  writer  by 

Professor  Brock. 

T      T 
Observe  that  starting  with  equation  (12),  CUU  C  =  W  and  pre- 

T 
multiplying  by  U  and  postmultiplying  by  U,  the  following  matrix  equa- 


tion is  obtained 


and 


T   T      T  T 
U  CUU  CU  =  U  W  U 
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T        T   T   3j 
U  CU  =  (U  W  UK 

T  T 
where  U  W  U  is  a  symmetric  matrix,  but  as  U  and  V  have  dimensions  (6xn) , 

where  n  is  the  number  of  pair  of  vectors  available,  this  formulation  of 

the  problem  requires  taking  the  square  root  to  an  (nxn)  symmetrix  matrix, 

with  the  disadvantage  of  dealing  with  a  big  matrix,  for  example  (1000x1000) 

Accordingly,  return  to  equation  (12) 

T      T 
CUU  C  =  W 

T 
Let  UU  =  A,  and  without  much  loss  of  generality  assume  that  the 

eigenvalues  of  A  are  distinct.   In  this  problem  A  is  a  (6x6)  matrix. 

Then 

A  =  I 

where  M  is  the  modal  matrix  of  A 


.V 


2 
and   L  is  the  spectral  matrix  of  A. 

2 
This  is  true  since  A  is  symmetrical,  L  is  simply  a  diagonal  matrix 

of  the  eigenvalues  of  A,  and  M  is  a  square  matrix  whose  columns  are  the 

T    -1 

corresponding  eigenvectors  of  A.  M  is  orthogonal  so  that  M  =  M 

2 
Since  L  is  a  diagonal  matrix,  it  is  easy  to  find  its  square 

root,  L.   Theoretically  there  are  a  number  of  different  choices  of  L 
since  either  the  positive  or  the  negative  square  root  of  each  eigen- 
values can  be  chosen,  but  in  this  application  it  is  convenient  and  satis- 
factory to  use  only  the  positive  values. 

Now  .   ,„2UT   „TT„T 

A  =  ML  M  =  MLLM 

T  T 
=  MLL  M 

T 

=  EE 

where  E  »  ML  is  (mxm) ,  i.e.,  6  x  6  in  the  applied  problem. 

Let  B  =  WT 
and  equation  (12)  now  takes  the  form 
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CEETC  =  B 

T   T      T 
E  CEE  CE  =  E  BE  =  D 

t 
where  D  is  an  (m^ra)  symmetric  matrix.   Thus, 

E  CE  =  D2 

Like  matrix  A,  D  can  be  written  in  terms  of  its  modal  and  spectral 

matrices 

2  T 
D  =  ml  m 

TT   =  mlm 


Therefore 

T       T 
E  CE  -  mlm 

Now  E"1  -  L"1!*"1  =  L'V" 

where  L  "  is  simply  the  diagonal  matrix  of  the  reciprocals  of  the  corres- 
ponding elements  of  L, 

Thus 

-1    T  -1  T 

C  ■  ML   mlm  L  M 

This  procedure  has  been  written  in  the  FORTRAN  IV  language  to  be 
used  in  the  IBM  System  360  computer,  and  in  the  actual  program  is  output 
3D. 
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SECTION  5 
SOLUTION  OF  EQUATION  FOR  CRITERION  4 
The  equation  (13),  that  minimize  the  criterion  4,  is  nonlinear  in 
the  unknown  matrix  C  (=  S  ) ,  and  an  iteration  method  is  used  for  its 
solution. 

Given  JB>      ,  choose  a  convenient  f  £.  &      .   Now  from  equation  (13), 
consider  the  iterative  scheme  indicated  by 

C.P  +  PC.  -  Q  =  f  S.(S.R  +  RS.  -  Q)S,  (17) 

where 

s.  -  c:\ 

i    i-I 
and      C  =  Zl  *    fS 

0        — 1  +J 

C    is  the  matrix  that  satisfies  criterion  1 

CL   is  the  matrix  that  satisfies  criterion  2 

Now,  having  C  ,  S  is  determined,  from  which  the  right  hand  side 

of  equation  (17)  can  be  evaluated,  and  the  following  equation  can  be 

obtained: 

C.P  +  PC.  =  0*  (18) 

where 

Q*  =  Q  +  f  S^S.R  +  RS.  -  Q)S. 
J   i   i      11 

Equation  (18)  is  of  the  same  form  as  equation  (10),  and  can  be 

solved  using  the  same  procedure. 

The  C,  obtained  from  equation  (18),  is  used  again  to  obtain  S .  . , 
i  1+1 

and  a  new  C.  .  is  obtained  using  the  iterative  procedure  until  two 
i+l 

successive  approximations  of  C.  become  less  than  a  specified  allowable 


error,  viz, 


«WJk  -  <Vjk    ■  A 
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where  ^   in  the  actual  program  is  taken  as  10~  ,  and  ()   implies 

J* 

comparing  element  by  element. 

When  the  convergence  is  achieved,  j*  is  increased  and  the  iterative 
procedure  is  carried  out  again.  This  procedure  is  followed  until  finally 
convergence  is  obtained  with  f  =  /3 

The  reason  for  using  values  of  j*  initially  less  than  A  is  to 
assure  and  optimize  convergence  of  the  entire  process.  For  example,  in 
one  study  where   /3  =0.8,  it  was  found  convenient  to  use  the  sequence 

f   -  0.1,0.2,0.3,. ..,0.8. 
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SECTION  6 
DIGITAL  COMPUTER  PROGRAMS 
Two  main  computer  programs  have  been  written  in  FORTRAN  IV  Language 
to  be  used  with  the  IBM  System  360  computer.   These  programs  are  called 
CECI  1  and  CECI  2.   They  make  use  of  several  subroutines  and  to  assure 
accuracy  the  main  programs  and  the  subroutines  employ  double  precision 
arithmetic. 

In  program  CECI  1,  the  following  subroutines  are  used: 
ADDEM  AINVM 

BIGCM  BIGSM 

BLANCA  JACVAT 

MAUSS3  OUTPU 

PRINM  RECIP 

REDEM  SQRERT 

SQROOT  SQUAM 

TRACM  TRANM 

ULTMA  ULTMB 

In  the  program  CECI  2,  the  following  subroutines  are  used: 
ADDEM  AINVM 

BIGCM  BIGSM 

BLANCA  MAUSS3 

OUTPUT  PRINM 

REDEM  SQUAM 

TRACM  TRANM 

ULTMA  ULTMB 

Both  programs  make  use  of  a  number  of  functions  from  the  monitor 
library  tape  which  are  not  listed  here. 
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Program  CECI  1  solves  the  equations  that  satisfy  criteria  numbers 

1,  2,  and  3.   The  program  gives  one  output  for  each  of  criteria  1  and 

2,  and  four  outputs  for  criterion  number  3.   For  each  of  these  outputs 
the  program  prints  the  compliance  matrix  C,  and  for  each  such  matrix, 
computes  and  prints  out  three  criteria  numbers,  as  follows: 

CRITERION  NUMBER  ONE     =ta  [(CU  -  V)T(CU  -  V)] 
CRITERION  NUMBER  TWO     =1jl  [(SV  -  U)T(SV  -  V)] 
CRITERION  NUMBER  THREE   »"frl  [(CU  -  V)T(U  -  SV) ] 
Program  CECI  2  solves  the  equation  that  satisfies  criterion  number 
4,  prints  the  compliance  matrix  C,  and  gives  four  criterion  numbers;  the 
first  three  are  the  same  as  in  CECI  1  and  the  fourth  is  given  by: 

CRITERION  NUMBER  FOUR    =t>lt(CU  -  V)T(CU  -  V)  +^8  (SV-U)T 

(SV-  U)] 

Most  of  the  matrices  are  regarded  as  doubly  subscripted  arrays  as 
follows: 

If  there  are  M  rows  and  N  columns,  the  array  in  the  computer  memory 
has  dimension  (M  +  1,  N)}  for  example:   if  the  matrix  has  6  rows  and  50 
columns,  the  dimension  of  the  array  is  (7  x  50).   This  extra  row  permits 
the  identification  of  the  matrix  for  debugging  purposes  and  gives  this 
information  to  those  subroutines  which  use  it  for  processing  or  monitor- 
ing purposes. 

Let  e.  .  refer  to  the  array  and  let  a.  .  refer  to  the  matrix.   Then 
i,J  i,J 

e    is  an  integer  identifying  the  array  and  the  matrix,   e    and  e.  - 
1,1  1,21,  j 

are  integers  M  and  N  respectively.   For  the  remainder  of  the  elements, 

e.  .,  .  ■  a.  ..   The  other  matrices  used  in  the  main  programs  that  have 

1+liJ    i»J 

the  dimensions  (N,N),  are  used  as  scratch  matrices.   They  are  needed 
when  calling  some  subroutines  that  have  dummy  subscripted  variables,  due 
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to  the  mode  of  storage  (column-wise)  used  by  the  IBM  System  360  computer. 
The  array  for  an  M-element  column  matrix  has  dimension  (M  +  1),  where 
the  first  element  is  the  integer  M. 

The  programs  can  be  used  in  two  different  ways:  (1)  for  classroom  or 
more  investigation  purposes  and  (2)  for  use  with  data  from  an  actual  ex- 
periment.  This  is  done  by  means  of  the  card  called  MODE,  where  MODE  -   1, 
2,  respectively. 

If  MODE  =  1,  the  program  will  read  the  matrices  C  and  U  from  the 
data  cards,  and  generate  matrix  V  =  CU.   For  testing  purposes  random 
errors  are  introduced  in  matrices  U  and  V.  The  explanation  of  how  these 
errors  are  introduced  will  be  given  in  the  next  section. 

If  MODE  =  2,  the  programs  will  read  the  matrices  U  and  V  from  the 
data  cards.   This  data  would  be  obtained  from  an  actual  experiment. 

The  explanation  of  how  to  use  these  programs,  and  how  the  data  is 
punched  in  cards  will  be  given  in  Appendix  A. 

Unless  other  FORMAT  statements  are  called  for,  the  only  printed  out- 

2 
put  is  as  follows: 

a)   In  main  program  CECI  1,  with  MODE  =  1. 

(1)  MATRIX  P  FOLLOWS 

MATRIX  NUMBER  1,  6  ROWS,  6  COLUMNS 
Listing  of  matrix  P  in  FORMAT  (IX,  8D12.4) 

(2)  MATRIX  BE  FOLLOWS 

MATRIX  NUMBER  2,  6  ROWS,  NN  COLUMNS 
Listing  of  matrix  BE  in  FORMAT  (IX,  8D12.4) 


2 
The  nomenclature  used  in  the  outputs  of  the  programs  CECI  1  and 
CECI  2  is  given  in  Appendix  C;  see  pages  76  and  89. 
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(3)  MATRIX  EB  FOLLOWS 

MATRIX  NUMBER  3,  6  ROWS,  NN  COLUMNS 
Listing  of  Matrix  EB  in  FORMAT  (IX,  8D12.4) 

(4)  ERROR  IN  BLANCA  =  (number  given) 

(5)  MATRIX  U  FOLLOWS 

MATRIX  NUMBER  4,  6  ROWS,  NN  COLUMNS 
Listing  of  matrix  U  in  FORMAT  (IX,  8D12.4) 

(6)  MATRIX  V  FOLLOWS 

MATRIX  NUMBER  5,  6  ROWS,  NN  COLUMNS 
Listing  of  matrix  V  in  FORMAT  (IX,  8D12.4) 

(7)  ELAPSED  TIME  =  (time  given  in  sec.) 

(8)  OUTPUT  NO.  1,  MINIMIZES  V  SQUARE  ERROR 
DESIRED  MATRIX  IS 

MATRIX  NUMBER  17,  6  ROWS,  6  COLUMNS 

Listing  of  desired  compliance  matrix  in  FORMAT  (IX,  8D12.4) 

CRITERION  NUMBER  1  EQUALS  (value  given) 

CRITERION  NUMBER  2  EQUALS  (value  given) 

CRITERION  NUMBER  3  EQUALS  (value  given) 

(9)  ELASPED  TIME  =  (time  given  in  sec.) 

(10)  OUTPUT  NO.  2,  MINIMIZES  U  SQUARE  ERROR 
DESIRED  MATRIX  IS 

MATRIX  NUMBER  33,  6  ROWS,   6  COLUMNS 

Listing  of  desired  compliance  matrix  in  FORMAT  (IX,  8D12.4) 

CRITERION  NUMBER  1  EQUALS  (value  given) 

CRITERION  NUMBER  2  EQUALS  (value  given) 

CRITERION  NUMBER  3  EQUALS  (value  given) 

(11)  ELAPSED  TIME  ■  (time  given  in  sec.) 
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(12)  OUTPUT  NO.  3A,  MINIMIZES  TRIANGLE  ERROR 
DESIRED  MATRIX  IS 

MATRIX  NUMBER  XXX,  6  ROWS,  6  COLUMNS 

(XXX  may  have  any  large  integer  value,  due  to  the  iteration 
procedure  used  to  obtain  output  3A) 
Listing  of  desired  compliance  matrix  in  FORMAT  (IX,  8D12.4) 
CRITERION  NUMBER  1  EQUALS  (value  given) 
CRITERION  NUMBER  2  EQUALS  (value  given) 
CRITERION  NUMBER  3  EQUALS  (value  given) 

(13)  ELASPED  TIME  -  (time  given  in  sec.) 

(14)  OUTPUT  NO.   3B,  MINIMIZES  TRIANGLE  ERROR 
Same  as  number  (12) 

(15)  ELAPSED  TIME  =  (time  given  in  sec.) 

(16)  OUTPUT  NO.  3C,  MINIMIZES  TRIANGLE  ERROR 
Same  as  number  (12) 

(17)  OUTPUT  NO.  3D,  MINIMIZES  TRIANGLE  ERROR 
Same  as  number  (12) 

(18)  ELAPSED  TIME  =  (time  given  in  sec.) 

(19)  TOTAL  ELAPSED  TIME  =  (time  given  in  sec.) 
b)  In  main  program  CECI  1,  with  MODE  =  2. 

(1)  MATRIX  U  FOLLOWS 

MATRIX  NUMBER  1,  6  ROWS,  NN  COLUMNS 
Listing  of  matrix  U  in  FORMAT  (IX,  8D12.4) 

(2)  MATRIX  V  FOLLOWS 

MATRIX  NUMBER  2,  6  ROWS,  NN  COLUMNS 
Listing  of  matrix  V  in  FORMAT  (IX,  8D12.4) 

(3)  ELAPSED  TIME  -  (time  given  in  sec.) 
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(4)   OUTPUT  NO.  1,  MINIMIZES  V  SQUARE  ERROR 
DESIRED  MATRIX  IS 

MATRIX  NUMBER  14,  6  ROWS,  6  COLUMNS 

Listing  of  desired  compliance  matrix  in  FORMAT  (IX,  8D12.4) 

CRITERION  NUMBER  1  EQUALS  (value  given) 

CRITERION  NUMBER  2  EQUALS  (value  given) 

CRITERION  NUMBER  3  EQUALS  (value  given) 

From  here  the  printed  out  is  the  same  as  when  MODE  =  1. 

c)   In  main  program  CECI  2,  with  MODE  =  1. 

(1)  MATRIX  P  FOLLOWS 

MATRIX  NUMBER  1,  6  ROWS,   6  COLUMNS 
Listing  of  matrix  P  in  FORMAT  (IX,  8D12.4) 

(2)  MATRIX  BE  FOLLOWS 

MATRIX  NUMBER  2,  6  ROWS,  NN  COLUMNS 
Listing  of  matrix  BE  in  FORMAT  (IX,  8D12.4) 

(3)  MATRIX  EB  FOLLOWS 

MATRIX  NUMBER  3,  6  ROWS,  NN  COLUMNS 
Listing  of  matrix  EB  in  FORMAT  (IX,  8D12.4) 

(4)  ERROR  IN  BLANCA  -  (value  given) 

(5)  MATRIX  U  FOLLOWS 

MATRIX  NUMBER  4,  6  ROWS,  NN  COLUMNS 
Listing  of  matrix  U  in  FORMAT  (IX,  8D12.4) 

(6)  MATRIX  V  FOLLOWS 

MATRIX  NUMBER  5,  6  ROWS,  NN  COLUMNS 
Listing  of  matrix  V  in  FORMAT  (IX,  8D12.4) 

(7)  ELAPSED  TIME  -  (time  given  in  sec.) 

(8)  ELAPSED  TIME  =  (time  given  in  sec.) 

(9)  NUMBER  OF  ITERATIONS  -  (number  given) 
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(10)  OUTPUT  NO.  4,  MINIMIZES  V  +  BETA  x  U  SQUARE  ERRORS 
DESIRED  MATRIX  IS 

MATRIX  NUMBER  XXX,  6  ROWS,  6  COLUMNS 

Listing  of  desired  compliance  matrix  in  FORMAT  (IX,  8D12.4) 

BETA  =  (number  given) 

CRITERION  NUMBER  1  EQUALS  (value  given) 

CRITERION  NUMBER  2  EQUALS  (value  given) 

CRITERION  NUMBER  3  EQAULS  (value  given) 

CRITERION  NUMBER  4  EQUALS  (value  given) 

(11)  ELAPSED  TIME  =  (time  given  in  sec.) 

(12)  TOTAL  ELAPSED  TIME  =  (time  given  in  sec.) 
d)  In  main  program  CECI  2,  with  MODE  ■  2. 

(1)  MATRIX  U  FOLLOWS 

MATRIX  NUMBER  1,  6  ROWS,  NN  COLUMNS 
Listing  of  matrix  U  in  FORMAT  (IX,  8D12.4) 

(2)  MATRIX  V  FOLLOWS 

MATRIX  NUMBER  2,  6  ROWS,  NN  COLUMNS 
Listing  of  matrix  V  in  FORMAT  (IX,  8D12.4) 
From  here  the  printed  out  is  the  same  as  when  MODE  =  1. 
Actual  printouts  for  both  programs,  when  MODE  =  1,  are  shown  in 
Appendix  C. 

The  listing  of  the  main  programs  and  subroutines  will  be  given  in 
Appendices  A  and  B  respectively. 
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SECTION  7 
TESTING  OF  THEORY  AND  DIGITAL  COMPUTER  PROGRAMS 
In  this  investigation  the  digital  computer  programs  written  to  solve 
the  matrix  equations  that  satisfy  the  four  criteria  of  optimality  have 
not  been  tested  using  actual  (mechanical)  experimental  data.   The  data 
have  been  generated  by  assuming  a  symmetric  positive  definite  matrix  C, 
and  a  set  of  vector  forces  [  F.>  in  the  form  of  matrix  U.   From  these  two 
matrices  the  corresponding  set  of  vector  deflections j  D.  |  in  the  form 
of  matrix  V  were  obtained  as  V  =  CU. 

To  simulate  experimental  data,  matrices  U*  and  V*  were  formed  in  the 
following  way: 

U*  =  U  +  E(U) 

V*  =  V  +  E(V) 

where  E(U)  and  E(V)  are  matrices  of  errors. 

Then  matrices  U*  and  V*  were  employed  as  if  they  were  composed  of 
experimental  data,  to  be  used  in  the  programs  to  recover  an  approximation 
matrix  C*  for  C. 

To  form  error  matrices  E(U),  and  E(V),  a  mathematical  procedure  was 
devised  to  perform  some  alterations  in  each  element  of  the  matrices,  U 
and  V,  in  a  random  manner  with  an  upper  limit  to  the  absolute  value  of 
the  percentage  of  error  being  introduced.  This  procedure  has  not  been 
tested  for  randomness,  but  the  method  used  to  obtain  the  alterations 
suggests  that  the  values  of  the  elements  in  matrices  E(U)  and  E(V)  are 
indeed  pretty  good  random  errors. 


3 
Note  that  the  use  of  the  asterisk,*,  is  different  in  this  section 

from  its  use  in  section  2. 
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The  random  errors  are  generated  by  a  subroutine  called  BLANCA. 
Because  the  first  row  of  most  of  the  matrices  are  used  for  identifying 
purposes,  the  subroutine  makes  use  of  the  element(l,l) ,  which  is  an 
identifying  integer  called  NUMBER. 

Subroutine  BLANCA  modifies  a  given  matrix  A  =  [  a.   ]  by  introduc- 
ing random  errors  E  =  [e  .],  resulting  in  a  modified  matrix  A*  =  A+E. 
One  specifies  the  maximum  percentage  error  to  be  applied. 

First  a  number  9    =  3.6327  +  (matrix  identification  NUMBER)  is 
formed.  ?    is  increased  for  each  element  by  adding  a  real  constant  (3.6327) 
It  is  not  allowed  to  be  greater  than  9.999;  this  is  done  by  subtracting  a 
real  constant  (8.4153)  each  time  f   becomes  greater  than  9.999. 

Next  an  exponent  6     ■  (I/10+J/25-K). 62832)  is  formed.   Next  the  num- 
ber  a  *  f         is  formed.   Next  the  number  *)  =  (fractional  part  of  d  )  is 
obtained.   Finally  the  relative  error  is  determined  as  to  x  r,  where  r  is 
the  preselected  maximum  relative  error.   Then 

e.  .  =  +  1)   x  r  x  a.  . 
1J   -  C        lj 

The  algebraic  sign  is  given  by  the  first  digit  in  the  decimal  ex- 
pression of  a  .   If  this  is  even,  the  sign  is  positive;  if  this  is  odd, 
the  sign  is  negative. 

A  special  treatment  is  used  if  a.  .  =  0,  for  which  e.  .  ■  +  7)     ,  the 
sign  being  determined  as  before.   This  has  obvious  disadvantages  if  the 
average  value  of  the  elements  in  A  is  small,  and  it  is  suggested  that 
future  work  make  an  appropriate  modification  to  the  present  procedure. 

In  retrospect  it  appears  that  perhaps  a  more  truly  random  variation 
could  have  been  attained  by  using  ~h    =  (fractional  part  of  100  a       ), 
say.   Furthermore,  it  would  have  been  desirable  to  incorporate  such  ran- 
domizing influences  as  the  date,  or  local  time,  or  job  sequence  which 
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should  become  available  internally  to  the  computer,  but  which  were  not 
available  at  the  time  this  work  was  done. 

For  convenience  in  comparing  values  of  C*  with  C,  to  see  how  good 
the  recovered  matrix  C*  was,  most  of  the  compliance  matrices  C  used  were 
diagonal  matrices.   However,  there  is  no  loss  in  generality  as  may  be 
seen  from  the  following  argument. 

Starting  with  the  matrix  equation  CU  =  V,  with  diagonal  C,  and  pre- 
multiplying  this  equation  by  an  orthogonal  matrix  M,  it  can  be  written  as 

MClAlU  =  MV  (19) 

T 
where  the  insertion  of  MM  in  equation  (19)  between  C  and  U  does  not  al- 

T 
ter  the  equality,  due  to  the  fact  that  M  M  =  I. 


Now  let 


MCMT  =  C 


MU  =  U 
MV  =  V 
equation  (19)  can  be  written  now  as 

CU  =  V  (20) 

Equation  (20)  has  the  same  form  as  original  matrix  equation  CU  =  V, 
but  with  the  difference  that  starting  with  C  (a  diagonal  matrix) ,  the 
new  symmetrical  matrix  C  is  no  longer  a  diagonal  matrix,  and  instead  all 
its  elements  are  generally  non-zero. 

Thus  a  general  problem  CU  =  V  could  be  written  in  the  particular 

form  CU  =  V,  with  diagonal  C,  if  only  the  orthogonal  matrix  M,  which  dia- 

T—       T-  T— 

gonalizes  C  were  known.  C  =  M  CM,  U  =  M  U,  and  V  =  M  V. 

The  particular  validity  of  this  argument  was  verified  by  two  types 

of  numerical  experiment.   In  the  first  type,  nondiagonal  C's  were  used, 

errors  were  introduced  in  U  and  V  ■  CU,  and  the  recovered  matrix  C*  was 
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successfully  compared  with  the  original  C.   In  the  second  type,  an  ortho- 
gonal matrix  M  was  generated,  and  matrices  C,  U,  and  V  were  formed  using 

M,  diagonal  C,  U,  and  V.   Then  random  errors  were  introduced  in  both  sets 

-  T 

of  data  (U,  V  and  U,V)  and  C*  and  C*  were  obtained.   Comparison  of  M  C*M 

with  C  was  quite  satisfactory. 

From  theory  the  minimum  number  of  independent  pairs  of  vectors  F  and 
D  needed  to  solve  the  problem  is  4.   This  is  because  in  the  compliance 
matrix  C  there  are  21  unknowns,  and  21  independent  equations  must  be  form- 
ed.  Thus  would  require  four  pairs  of  6-vectors  F  and  D  since  the  notion 
of  3%  pairs  is  excluded.   However,  the  particular  methods  employed  in  the 
programs  seem  to  demand  a  larger  number.   The  minimum  number  of  pairs  of 
vectors  needed  to  obtain  "good"  results  was  six. 

In  one  case,  a  problem  with  24  pairs  of  vectors  was  solved  twenty 
times,  using  a  maximum  random  error  of  ten  percent.   The  errors  were  in- 
troduced in  a  different  manner  for  each  of  the  solutions.   From  the  twenty 
compliance  matrices  obtained  in  this  way,  an  average  compliance  matrix 
was  obtained.   This  matrix  looks  "very  good",  when  compared  element  by 
element  with  the  original  compliance  matrix.   This  indicates  that  if  the 
value  of  n  is  so  large  as  to  overload  computer  memory  or  as  to  prove  in- 
convenient for  loading  data,  the  U  and  V  matrices  may  be  conveniently 
partitioned,  a  number  of  "smaller  problems"  dealt  with,  and  the  answers 
averaged.   One  could  regard  the  case  described  as  forming  a  single  pro- 
blem with  n  =  480  which  was  actually  treated  as  the  average  of  the  re- 
sults of  twenty  problems,  each  having  n  =  24. 
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SECTION  8 
CONCLUSIONS  AND  RECOMMENDATIONS 

There  does  not  seem  to  be  any  direct  way  of  testing  the  validity 
of  the  theory  and  the  computer  programs  which  implement  it  other  than  by 
taking  actual  examples.   These  examples  were  constructed  by  assuming 
values  for  matrices  C  and  U  and  generating  matrix  V  =  CU.  Then,  in  most 
cases,  errors  were  applied  so  as  to  get  modified  matrices  U*  and  V*. 
Then,  from  these  matrices,  considered  as  input  data  to  the  program,  the 
various  procedures  were  used  to  recover  an  optimal  value,  C*,  which  is 
then  compared  with  the  original  C.   No  useful  or  convenient  method  has 
been  found  for  assessing  the  accuracy  of  the  recovery  other  than  by  com- 
paring term  by  term,  and,  to  simplify  this  task,  in  most  cases  C  was 
assumed  to  be  a  diagonal  matrix.   It  was  shown  in  Section  7  that  there  is 
no  loss  in  generality  in  doing  this. 

One  of  the  things  that  might  be  expected  to  influence  the  ease  and 
accuracy  of  the  recovery  is  the  relation  between  the  eigenvectors  of  the 
matrix  C.   Tests  were  conducted,  at  the  one  extreme,  with  diagonal  C 
having  identical  diagonal  elements,  and  at  the  other  extreme,  with  dia- 
gonal matrices  having  diagonal  elements  differing  by  a  factor  of  3125. 
In  all,  fairly  extensive  numerical  work  was  investigated  and  various 
maximum  percentage  errors  were  introduced.   The  output  produced  as  a 
result  of  these  studies  is  much  too  volumnious  to  be  reported  in  detail 
and  it  has  not  been  found  possible  to  reduce  the  volume  of  information 
in  a  satisfactory  quantitative  way.   Accordingly,  it  is  necessary  to 
restrict  remarks  to  those  of  a  general  nature  which  will  be  made  in  the 
following.   In  general,  it  may  be  concluded  that  the  theory  and  the  pro- 
gram have  been  proved  to  be  accurate  and  effective,  and  in  combination 
they  provide  a  useful  tool  for  analysing  the  data  from  certain  types  of 
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structural  and  mechanical  experiments. 

The  following  matrices  have  been  used  as  compliance  matrices  C 
from  which  the  matrices  C-    v/ere  recovered. 
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6-_ 

0    0    0  3125^ 
Gives  the  desired  compliance  matrix  that  satisfies 
criterion  one.  With  this  procedure  the  types  of 
matrices  numbered  1,  4,  5,  and  6,  have  been  recovered 
quite  successfully  using  up  to  35%  maximum  random 
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errors.  Matrices  numbered  3  and  7  were  recovered  with 
up  to  5%  maximum  random  errors,  and  matrix  2  was 
recovered  with  errors  up  to  25%.  (Where  no  errors 
were  introduced,  each  original  was  accurately  recover- 
ed). 
OUTPUT  No.  2        Gives  the  desired  compliance  matrix  that  satisfies 

criterion  two.   All  the  matrices  were  recovered  same 
as  the  originals  when  no  errors  were  involved.  Matri- 
ces numbered  1,  2,  3,  4,  and  5  have  been  recovered 
quite  successfully  using  up  to  35%  random  errors, 
and  matrices  numbered  6  and  7  with  up  to  25%  random 
errors . 

OUTPUT  Nos.  3A,3B,  Gives  the  desired  compliance  matrix  that  satisfies 
3C 

criterion  three,  using  the  procedure  that  obtains  the 

square  root  of  a  non-symmetric  matrix.   All  the  matri- 
ces were  recovered  same  as  the  originals  when  no  errors 
were  involved.  With  matrices  numbered  1  to  6  this 
procedure  has  been  worked  successfully  with  up  to  357o 
random  errors,  and  with  matrix  numbered  7  with  up  to 
207o  random  error. 
OUTPUT  No.  3D       Gives  the  desired  compliance  matrix  that  satisfies 

criterion  three,  using  the  procedure  that  takes  square 
roots  of  two  different  symmetric  matrices.   Matrix 
numbered  1  is  the  only  one  that  has  not  been  recovered 
same  as  the  original  when  no  errors  were  involved. 
In  the  light  of  other  results  obtained,  this  failure 
to  recover  matrix  1  to  a  very  high  degree  of  accuracy 
is  unexpected  and  it  has  not  been  satisfactorily 
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explained.   Matrices  numbered  1  to  6  using  up  to  35% 
random  error,  have  been  recovered  quite  successfully, 
and  matrix  numbered  7  with  up  to  257o. 
OUTPUT  No.  4       Gives  the  desired  compliance  matrix  that  satisfies 

criterion  four.   Using  any  value  for  BETA,  when  no 
errors  were  involved,  all  the  matrices  were  recover- 
ed same  as  the  originals.  When  the  error  was  increased, 
better  results  were  obtained  for  small  BETA's.  The 
same  remarks  given  for  OUTPUT  No.  1,  applies  to  this 
output. 
To  have  a  better  idea  how  these  results  look,  see  Appendix  C,  where 
the  computer  solutions  for  a  particular  problem  with  zero  and  ten  percent 
random  errors,  are  given. 

The  work  presented  here  is  just  the  beginning  and  not  the  last  word 
in  this  new  area;  more  work  and  investigation  should  be  done.  Some  sug- 
gestions are  given  below: 

1)  Solution  of  equations  that  satisfy  criteria  1,  2,  and  4, 
can  be  investigated  using  the  procedures  given  by  Bickley  and  McNamee 
[4],  or  possible  modifications  to  take  advantage  of  the  fact  that  matrices 
A  and  B  are  equal. 

2)  Investigate  how  large  maximum  random  error  percentages  one 
can  give  to  matrices  U  and  V.  Give  to  these  matrices  different  maximum 
percentages  of  errors,  e.g.,  to  matrix  U,  5%,  and  to  matrix  V,  15%. 

3)  There  exists  a  fifth  criterion  given  by  Brock  [1],  that  has 
not  been  investigated  yet. 

4)  Devise  a  new  iterative  procedure  for  the  equation  that 
satisfies  criterion  four. 
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5)  Devise  a  convenient  criterion  which  measures  the  "overall" 
accuracy  of  recovery  of  a  matrix  C.  At  present,  one  merely  compares 
element  by  element. 

6)  Employ  advanced  and  sophisticated  statistical  theory  for  the 
analysis  of  errors  and  the  accuracy  of  recovery.   This  would  seem  to 
imply  constructing  an  adequate  set  of  definition  and  the  use  of  probabil- 
ity spaces,  etc.   The  present  analysis,  in  which  general  conclusions  are 
attempted  from  specific  cases,  leaves  much  to  be  desired. 

As  things  stand  now,  the  engineering  analyst  has,  in  this  theory 
and  program  implementation,  a  four  edged  tool  which  can  be  used  for  the 
analysis  of  experimental  data.   The  decision  of  which  criterion  yields 
the  "best"  result  is  one  which  must  be  made  on  the  basis  of  experience 
with  the  physical  problem. 


42 


BIBLIOGRAPHY 

1.  Brock,  John  E.   Second  Progress  Report,  contract  N161-26327  with 
John  F.  Wax,  Inc.,  and  United  States  Naval  Academy,  May  15,  1966. 

2.  Martin,  Harold  C.   Introduction  to  Matrix  Methods  of  Structural 
Analysis,  McGraw-Hill,  1966. 

3.  Frazer,  R.  A.,  Duncan,  W.  J.,  and  Collar,  A.  R.   Elementary  Matrices, 
Cambridge  University  Press  and  the  Macmillan  Co. ,  Cambridge  and  New 
York,  1947. 

4.  Bickley,  W.  G. ,  and  McNamee,  J.  Matrix  and  other  direct  methods  for 
the  solution  of  systems  of  linear  difference  equations,  Philosophi- 
cal Trans,  of  the  Royal  Society  of  London,  Series  A.,  No.   1005,  v. 
252,  21  January  1960:  69-131. 

5.  Larson,  H.  J.   Least  squares  estimation  of  the  components  of  a  sym- 
metrical matrix.   Techonometrics ,  v.  8,  No.  2,  May  1966:  360-362. 

6.  Brock,  John  E.   Third  Progress  Report,  contract  N161-26327  with 
John  F.  Wax,  Inc.,  and  United  States  Naval  Academy,  July  15,  1966. 

7.  Pulay,  Peter.   An  Iterative  Method  for  the  Determination  of  the 
Square  Root  of  a  Positive  Definite  Matrix.   Zeit.  d.Ang.Math.  u. 
Mech.  (ZAMM),  v.  46,  1966:  151. 


43 


APPENDIX  A 
INSTRUCTIONS  FOR  USE  AND  LISTING  OF  PROGRAMS 
A-l  Preparation  of  Data  Deck  and  Final  Assembly  of  Programs. 

The  preparation  of  the  data  deck  is  the  same  for  programs  CECI  1 
and  CECI  2,  but  has  some  alterations  depending  on  whether  MODE  =  1  or 
MODE  -  2. 

In  what  follows  the  instructions  of  how  to  prepare  the  data  deck  is 
given  for  both  types  of  modes. 

The  first  card  is  bbbbbbb6NNbbblbb ,  where  b  =  blank,  the  6  is 

the  number  of  rows,  NN  is  the  number  of  columns,  and  1  means  first  matrix 
to  be  read.   When  MODE  =  1,  NN  ■  06,  is  the  number  of  columns  in  matrix 
C;  when  MODE  =  2,  NN  varies  from  06  to  50,  and  is  the  number  of  columns 
in  matrix  U.   (Thus  U  is  limited  to  a  maximum  of  50  columns). 

Following  this  card,  comes  the  elements  of  matrix  C  or  U  as  the 
case  may  be,  punched  in  five  entries  to  a  card,  in  FORMAT  (5F14.6). 
Entries  for  the  first  row  are  thus  entered  until  all  elements  for  the  row 
have  been  punched.   Then,  starting  with  a  new  card,  entries  for  the  sec- 
ond row  are  punched,  and  so  on,  until  the  entries  for  the  sixth  row  have 
been  entered.   Each  row  starts  with  a  new  card. 

Card  bbbbbbb6NNbbb2bb ,  follows.   Then  the  entries  for  the 

matrix  U  (if  MODE  =  1)  or  matrix  V  (if  MODE  =  2)  are  punched  in  the  same 
manner  as  before. 

Example:   Punch  cards  as  in  Figure  A-l,  for  the  following  6x6  diag- 
onal matrix  as  C  matrix. 


C  - 


1. 
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6 
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Sixth  row 

Fifth  row 

Fourth  row 

Third  row 

Second  row 
First  row 

First  card 


NOTE:  b  -  blank,  B  -  14b 's. 


FIGURE  A-l 

Note  that  element  C(l,l)  is  punched  between  fields  1  to  14,  element 
C(2,2)  between  fields  15  to  28,  and  so  on.   Note  too  that  each  row  has 
two  cards  because  there  are  six  elements  for  a  row,  and  in  each  card 
only  five  consecutive  elements  can  be  punched. 

If  more  than  one  set  of  data  is  used,  the  second  set  comes  after 
the  last  card  of  the  first  set,  starting  again  with  the  first  card 

bbbbbbbNNbbblbb 

and  so  on. 

In  order  to  stop  the  program,  the  final  card  of  the  data  deck  is 
punched  as  follows  bbbbbbb0b6bbb6bb 

The  complete  program  deck  includes?  control  cards,  main  program, 
subroutine  package,  and  data  deck.   The  assembly  of  these  components  is 


45 


shown  in  Figure  A-2,  where  the  first  control  card  is  from  the  CECI  1 
program.   To  use  program  CECI  2,  change  the  word  CECI  1  to  CECI  2. 


CONTRDL\  /* 
CARDS  \ 


FIGURE  A-2 
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A-2  Listing  of  Program  CSCI  1. 

The  cards  of  program  C^CI  1  are  numbered  from  0000  to  0159. 
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A-3  Listing  of  Program  CBCI  2. 

The  cards  of  program  CECI  2  are  numbered  from  00000  to  01570. 
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APPENDIX  B 
SUBROUTINE  PROGRAMS 
B-l  Description  of  Subroutines. 

SUBROUTINE  ADDEM  (A,B,C,D,M,L)  provides  the  result   of  A  =  BfDC, 
where  D  is  a  scalar  constant,  and  A,  B,  and  C  are  double  precision  matri- 
ces dimensioned  (M,L) ,  where  M  and  L  are  not  less  than  3.  With  D  =  1.0, 
A  -  B4C;  with  D  -  -1.0,  A  =  B  -  C. 

The  matrices  B  and  C  from  the  calling  program  must  have  identical 
order,  otherwise  a  message  is  printed  and  the  execution  is  halted. 

SUBROUTINE  AINVM  (A,B,C,D,L,N)  forms  the  inverse  A  =  B"!    of  input 
matrix  B.  Matrices  A,  B,  C,  and  D  are  double  precision,  with  A  and  B 
dimensioned  (L,N).  Matrices  C  and  D,  which  are  used  as  scratch  matrices, 
have  dimension  (N,N),  where  L  =  N+l,  and  N  is  not  less  than  3. 

A  message  is  printed  if  the  input  matrix  B  is  non-square  and  the 
execution  is  halted. 

SUBROUTINE  BIGCM  (C0L,Q)  forms  the  21-element  vector  COL  from  the 
6x6  matric  Q,  according  to  the  scheme  described  in  Section  3;  if  the 
input  matrix  Q  is  non-square  a  message  is  printed  and  the  execution  is 
halted. 

SUBROUTINE  BIGSM  (BIG,P)  forms  the  21  x  21  matrix  BIG,  from  the  6x6 
matrix  P  according  to  the  scheme  described  in  Section  3;  if  the  input 
matrix  P  is  non-square  a  message  is  printed  and  the  execution  is  halted. 

SUBROUTINE  BLANCA  (C0,BE,N)  forms  matrix  CO  from  matrix  BE  by  intro- 
ducing random  errors  in  the  latter  according  to  the  method  described  in 
Section  7 .   The  arrays  CO  and  BE  are  dimensioned  (7 ,N) .   The  maximum 
relative  error,  denoted  by  the  symbol  r  in  the  exposition  in  Section  7, 
is  governed  by  the  computer  variable  PERC,  which  is  one  tenth  the  relative 
error.   If  PERC  =  l.D-2  -  .0100...,  then  r  =  .100...,  and  a  maximum  of 
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10%  error  is  introduced.   PERC  may  be  changed  by  changing  card  number 
0313  in  the  subroutine. 

SUBROUTINE  JACVAT  (B,N,N0YES,  RR,SS,EIRU,NDIM)  gives  the  eigenvalues 

and  eigenvectors  of  the  symmetric  matrix  B.   If  B  is  non- symmetrical  the 

T 
subroutine  first  makes  B  symmetric  by  forming  %B+%B  .   RR  is  the  spectral 

matrix  that  contains  the  eigenvalues  of  B.   SS  is  the  modal  matrix  of  B 

that  contains  the  eigenvectors.   EIRU  is  the  square  root  of  the  elements 

in  the  diagonal  matrix  RR.   N  is  the  dimension  of  B.   NOYES  is  0  or  1; 

if  0,  no  eigenvalues  will  be  recovered;  if  1,  the  eigenvalues  will  be 

recovered  in  matrix  RR.   In  the  application  made  in  this  thesis,  NDIM  is 

the  same  as  N.   All  the  matrices  and  variables  are  in  double  precision. 

SUBROUTINE  MAUSS3  (N,EP,A,X,KER)  is  a  double  precision  FORTRAN  IV 
language  subroutine  which  inverts  a  symmetric  matrix.  A  is  the  matrix 
to  be  inverted,  X  is  its  inverse.   Both  have  dimension  of  (N,N).   EP  is 
a  tolerable  error,  taken  as  l.D-20  in  the  present  work.   KER  is  an  output 
flag  having  the  value  1,  if  all  has  gone  well,  or  2,  if  there  has  been 
any  difficulty  of  execution. 

SUBROUTINE  OUTPU  (R,B,U,V,T1,T2,T3,T5,LL,L)  has  inputs  R,  B,  U,  V, 
where  R  is  the  desired  compliance  matrix  .   B  is  the  desired  stiffness 
matrix,  and  U,  V  are  the  data  for  the  problem  (such  that  RU  =  V,  approxi- 
mately).  It  prints  the  matrix  R,  and  computes  and  prints  the  three 
criterion  numbers  o(^  (i  =  1,2,3).   Matrices  T1-T5  that  appear  in  the 
calling  statement  are  scratch  matrices.  LL  and  L  are  dummy  variables 
used  in  the  dimension  of  the  matrices  and  are  related  to  the  number  n  of 
vector  pairs  used. 

SUBROUTINE  PRINM  (A,L)  causes  the  (6,L)  matrix  A  to  be  printed, 
where  L  is  not  less  than  3.   The  first  line  printed  contains  the  fol- 
lowing information:   MATRIX  NUMBER  XX,  YY  ROWS,  ZZ  COLUMNS.   This  is 
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followed  by  a  print  of  the  matrix,  line  by  line,  in  FORMAT  (IX,  8D12.4). 

SUBROUTINE  RECIP  (A,B)  A  and  B  are  diagonal  matrices  (6,6)  and  A  = 
B"  ,  where  B  is  the  input  matrix. 

SUBROUTINE  REDEM  (A,IDENT)  reads  into  the  computer  memory  from  the 
data  deck  which  is  punched  as  described  in  Appendix  A.  A  is  the  "name" 
of  the  matrix  and  IDENT  is  an  identifying  integer  to  assure  that  the 
card-group  which  is  read  is  indeed  the  one  desired. 

SUBROUTINE  SQRERT  (A,B,C)  solves  the  matrix  equation  that  satisfies 

criterion  3,  using  the  procedure  that  twice  takes  the  square  root  of  a 

T       T 
symmetric  matrix.   B  =  UU  ,  C  =  W  are  the  inputs,  and  A  is  the  desired 

compliance  matrix. 

SUBROUTINE  SQROOT  (A,CC,G)  takes  the  square  root  of  a  non-symmetric 
matrix  by  iteration  procedure.  Matrices  CC  and  G  are  inputs.   CC  is  the 
approximate  square  root;  G  is  the  matrix  whose  square  root  is  desired, 
and  A  is  the  desired  square  root. 

SUBROUTINE  SQUAM  (A,B)  forms  a  6  by  6  matrix  A  from  the  21-element 
column  vector  B,  arranging  the  results  properly  as  implied  in  Section  3. 

SUBROUTINE  TRACM  (T,A)  forms  the  trace  T  of  the  6  by  6  matrix  A. 

SUBROUTINE  TRANM  (A,B,LL,L)  forms  the  transpose  of  matrix  B(7,L), 
obtaining  matrix  A(LL,6),  where  L  is  no  less  than  3,  and  LL  =  L  +  1. 

SUBROUTINE  ULTMA  (A,B,C,MM,M,L,LL)  forms  the  matrix  product  A  =  BC. 
A(MM,M)  =  B(MM,L)  x  C(LL,M),  where  MM,  M,L,LL,  can  be  any  numbers  not  less 
than  3  for  which  the  matrices  are  conformable  for  multiplication.   A 
message  is  printed  if  the  matrices  B  and  C  are  not  conformable  and  execu- 
tion is  halted. 

SUBROUTINE  ULTMB  (A,B,C,LL,L)  multiplies  the  square  matrix  B  by  a 
column  matrix  C,  to  obtain  the  column  matrix  A.  A(LL)  =  B(LL,L)  x  C(LL), 
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where  L  is  any  number  not  less  than  3,  and  LL  =  L  +  1. 

SUBROUTINE  OUTPUT  (R,B,U,V,T1,T2,T3,T5,LL,L,BETA)  is  the  same  as 
SUBROUTINE  OUTPU.  The  only  differences  are  the  inclusion  of  parameter 
BETA  in  the  calling  statement,  and  the  printing  of  four  criterion  numbers 
instead  of  three,  o^-L   (i  =  1,2,3,4). 


B-2  Subroutine  Listing 

SUBROUTINE 

ADDEM 

AINVM 

BIGCM 

BIGSM 

BLANCA 

JACVAT 

MAUSS3 

OUTPU 

PRINM 

RECIP 

REDEM 

SQRERT 

SQROOT 

SQUAM 

TRACM 

TRANM 

ULTMA 

ULTMB 

OUTPUT 


CARDS 

NUMBERED 

FROM 

TO 

0160 

0182 

0183 

0208 

0209 

0233 

0234 

0288 

0289 

0336 

0337 

0481 

0482 

0532 

0533 

0559 

0560 

0576 

0577 

0594 

0595 

0616 

0617 

0638 

0639 

0692 

0693 

0727 

0728 

0744 

0745 

0763 

0764 

0791 

0792 

0814 

20000 

20360 
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APPENDIX  C. 

SAMPLE  PROBLEM  SOLVED  BY  COMPUTER 

C-l  Sample  Problem  Using  Program  CECI  1  with  0%   and  107o  Maximum  Random 
Errors . 

The  nomenclature  used  in  the  output  of  the  program  is  as  follows: 


MODE  =  1 
Matrix  P 
Matrix  BE 
Matrix  EB 
Matrix  U 
Matrix  V 
Output  No.  1 

Output  No.  2 

Output  No.  3A 


Output  No.  3B 


Output  No.  3C 


Output  No.  3D 


Data  will  be  generated  in  the  program. 

Original  matrix  C  as  read  from  cards. 

Matrix  U  as  read  from  cards. 

Matrix  V  from  V  =  CU. 

Matrix  BE  with  random  errors  introduced. 

Matrix  EB  with  random  errors  introduced. 

The  desired  compliance  matrix  which  satisfies 
criterion  one. 

The  desired  compliance  matrix  which  satisfies 
criterion  two. 

The  desired  compliance  matrix  which  satisfies 
criterion  three,  using  the  iterative  procedure  for 
taking  the  square  root  of  the  non-symmetric  matrix 
WTUUT. 

The  desired  compliance  matrix  which  satisfies 
criterion  three,  using  the  iterative  procedure  of 
taking  the  square  root  of  the  non-symmetric  matrix 

uuTwT. 

The  desired  compliance  matrix  which  satisfies 
criterion  three.   It  is  the  average  of  outputs 
3A  and  3B. 

The  desired  compliance  matrix  which  satisfies 
criterion  three,  using  the  procedure  that  twice 
takes  square  roots  of  symmetric  matrices. 


The  printed  output  for  the  problem  with  zero  error  appears  on  pages 
77  to   82  inclusive;  that  for  the  problem  with  10%  maximum  error  ap- 
pears on  pages  83   to  88   inclusive. 
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C-2  Sample  Problem  Using  Program  CECI  2  with  0%  and  10%  Maximum  Random 
Errors. 

The  nomenclature  used  in  the  output  of  the  program  is  as  follows: 


M0DE*1 
BETA  0.8 
Matrix  P 
Matrix  BE 
Matrix  EB 
Matrix  U 
Matrix  V 


Data  will  be  generated  in  the  program 

(The  sequence  used  is  P=  .1,.2,...,.8.) 

Original  matrix  C  as  read,  from  cards. 

Matrix  U  as  read  from  cards. 

Matrix  V  from  V  ■  CU. 

Matrix  BE  with  random  errors  introduced. 

Matrix  EB  with  random  errors  introduced. 


Output  No.  4   The  desired  compliance  matrix  which  satisfies 
criterion  four. 

The  printed  output  for  the  problem  with  zero  error  appears  on 

pages  90  to   93  inclusive;  that  for  the  problem  with  10%  maximum 

error  appears  on  pages  94  to  97  inclusive. 
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APPENDIX  1 
BRIEF  TREATMENT  OF  RELATED  PROBLEMS 
D-l  Case  whjre  there  are  no  Restrictions  on  C. 

In  the  development  of  this  thesis,  starting  with  the  matrix  equa- 
tion CU  =  V,  matrix  C  was  restricted  to  be  symmetric  and  positive  definite. 
Suppose  now  that  these  restrictions  should  be  removed,  and,  indeed,  let 
it  only  be  required  that  C  have  the  correct  dimensions  to  relate  matrices 
U  and  V  which  may  have  a  different  number  of  rows. 
Let  us  write  equation  (4)  in  a  modified  form 

oC  =  Dl   [(CU-V)T£(CU-V)]  (D-l) 

where  Z  is  arbitrary  and  might  be  thought  of  as  representing  a  set  of 
weighting  factors.   If  Z  is  the  unit  matrix,  equation  (D-l)  reduces  to 
equation  (4) . 

From  equation    (D-l) 

o<  -  ti  [UTCTZ  CU  -UTCT2V-VT2CU  +YTZV] 

T 
Note  that,  in  general,  C  zfc    C  since  the  requirement  of  symmetry 

is  not  made  here. 

Using  *  to  denote  partial  differentiation  with  respect  to  any  element 

x,  as  before,  one  gets 

<x*=ftL[uV*£CU  +uTcTzc*u  -uTcT*zv-\/T£C*u  ] 

-~1h   [  (CUUT-VUT)  (CT*J  (*  +*T  )]  (D-2) 
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Thus  o(,  is  made  stationary  with  respect  to  -^ach  element  of  C  if  C 

satisfies  the  equation 

T     T 
CUU  =  VU 

If  U  is  of  rank  m,  whera  mis  the  number  of  rows  in  U,  then 

c  -  vu^uirV1 

It  is  interesting  to  note,  that  the  result  is  independent  of  the 
arbitrary  weighting  matrix  Z. 

D-2  Case  where  C  is  Orthogonal. 

This  problem  arises  in  the  treatment  of  the  orientation  of  a  rigid 
body. 

If  C  is  orthogonal,  then 

T    -1 
C  -  C   =  Y 

where  the  notation  Y  =  C   is  used  for  convenience  in  manipulation. 

As  before  let  x  be  any  element  of  C,  and  let  *  denote  -~- 


Then  the  relation  CY  ■  I,  leads  to  C*  =  -CY*C,  and  equation  (4)  can  be 

written  as 

OC-#L  (UTU-U"rYV-VTCU    +V"V) 

so  that 

<**=&.  (-UTY*V  +  V^CY^CU) 

=  "^[(Y*)(CUVTC-VU"r)]  (D-3) 

Thus,  o(  is  stationary  with  respect  to  variation  in  any  of  the  ele- 
ments of  C,  if  C  satisfies  the  equation 

T      T 
CUV  C  =  VU 

which  will  be  so  if  C  is  any  of  the  several  values 

C  =  (vuTuvT)^(uvT)"1 
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T  T 
VU  UV  is  symmetric  matrix,  and  its  square  root  can  be  obtained  as 

described  in  Section  4. 

These  results  were  obtained  by  Professor  Brock  and  are  given  here 

so  as  to  collect  in  one  place  information  on  this  general  matrix  problem. 
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